arXiv:nucl-th/0211083v2 4 Dec 2002 


1 


Important experimental observables at RHIC * 

D. Kharzeev’^'l’ 

’^Nuclear Theory Group, 

Physics Department, 

Brookhaven National Laboratory, 

Upton, New York 11973-5000, USA 

In this talk I discuss the signihcance of the hrst RHIC measurements for establishing 
the properties of hot and dense QCD matter and for understanding the dynamics of the 
theory at the high parton density, strong color held frontier. Hopes and expectations for 
the future are discussed as well. 

1. Introduction 

RHIC began operation in 2000, culminating over ten years of development and con¬ 
struction and a much longer period of theoretical speculations about the properties of 
hot QCD matter produced in nuclear collisions in the collider regime. This Conference 
presents a good occasion to look at the hrst experimental results and to discuss their 
meaning and signihcance. It is also time to think about the important questions which 
at present remain unanswered, and the future studies necessary to answer them. 

RHIC is a machine entirely dedicated to the study of Quantum Chromo-Dynamics - 
the theory of strong interactions. Asymptotic freedom of QCD [i. [I ensures that 
the dynamics of quarks and gluons at sufficiently high density can be addressed by weak 
coupling methods. This includes both the thermalized quark-gluon plasma at high tem¬ 
perature and the wave functions of the colliding nuclei described, at small Bjorken x, by 
parton saturation and the Color Glass Condensate What have we learned 

about these systems from experiment? Has the dense quark-gluon matter been produced 
at RHIC? I will return to these questions at the end of the talk; now, let us discuss what 
the experiment is telling us. 

2. Looking at the first RHIC data 
2.1. Global observables 

Global observables are the most general characteristics of the collision, including par¬ 
ticle multiplicity, its dependence on the centrality of the collision and on rapidity, and 
azimuthal distribution of the produced particles. The centrality is determined by the 
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impact parameter in the collision; since this is not a quantity which can be measured 
directly, centrality is usually determined in terms of a certain cut in the multiplicity dis¬ 
tribution; e.g., a 0 — 10% centrality cut means that out of a given sample, 10% of the 
events which have the highest multiplicity have been selected. 

It is convenient to characterize centrality 
in terms of the number of “participants” - 
the nucleons which underwent an inelastic 
interaction in a given collision. Glauber 
theory [ can be used to correlate a cer¬ 
tain centrality cut with an average number 
of participants (for an explicit set of formu¬ 
lae for nuclear collisions, see e.g. [i. [0). 

This procedure can be independently ver- 
ihed by measuring the energy carried for¬ 
ward by spectator neutrons; for that pur¬ 
pose all RHIC experiments are equipped by 
Zero Degree Calorimeters. 
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Figure 1: Centrality dependence of the charged 
particle multiplicity near mid-rapidity in Au -I- Au 
collisions at = 20,130,and 200 GeV, from [ |l7| . 


2.1.1. Multiplicity 

Multiplicity in heavy ion collision tells 
us which fraction of the collision energy is 
inelastically transferred to secondary par¬ 
ticles. 

Theoretical expectations for hadron mul¬ 
tiplicities at RHIC varied by factor of five 
(see, e.g., [0), and the experimental ver¬ 
dict was thus eagerly awaited. After the 
commissioning of the machine, the hrst 
multiplicity results did not take long to 
come; they are shown on Fig. Q in com¬ 
parison to the multiplicity previously mea¬ 
sured in pp and pp collisions. 

The measured multiplicity appeared 
much smaller than most theoretical pre¬ 
dictions. Is this disappointing? To an¬ 
swer this question, let us recall that, by 
the very definition, an incoherent superpo¬ 
sition of independent nucleon-nucleon col¬ 
lisions yields multiplicity equal to the number of collisions times the multiplicity in NN 
collision. This trivial statement holds also in the presence of elastic rescatterings. Indeed, 
according to so-called AGK cutting rules 
cross section is given by 
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Figure 2: Centrality dependence of the charged 
particle multiplicity near mid-rapidity in Au -I- Au 
collisions at -/s = 130 GeV and 200 GeV, from [p^. 
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where the nuclear overlap function is 
TAB{b)= f d^sTA{s)TB{b-s), 


( 2 ) 


and the nuclear thickness function T/i(6) = dzpA{b, z) is the integral over the nuclear 
density. Integration over impact parameter b yields 


and correspondingly the particle multiplicity would scale as 
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Using the numbers of collisions (~ 1050) and participants (~ 340) in central (0 — 6% 
centrality cut) An An collisions from Glauber model calculations [0, we would thus 
conclude that AuAu multiplicity per participant pair should exceed NN multiplicity by 
factor of 3. Instead, the data at the highest RHIC energy of ^/s = 200 GeV show the 
difference of only about 50%. Given that any inelastic rescatterings in the hnal state 
can only increase the multiplicity, we therefore have an experimental proof of a high 
degree of coherence in multi-particle production in nuclear collisions at RHIG energies. 
The diagrams which allow to evade the AGK theorem are Gribov’s “inelastic shadowing” 
corrections [ 0| which correspond to the excitation of high-mass intermediate states in 
multiple scattering; the process thus no longer can be decomposed in terms of elementary 
NN interactions. In parton language, these contributions correspond to multi-parton 
coherent interactions. 


2.1.2. Centrality dependence 

The dependence of multiplicity upon the number of participants discussed above can 
be established by selecting different centrality cuts. The result is shown in Fig.§ one can 
see that the multiplicity per participant pair increases with centrality, but not quite as 
fast as it would if the NN collisions were independent. 

If we decompose the multiplicity measured in A^A^ collisions at some energy ^/s into 
a fraction X(s) coming from “hard” processes, and the remaining fraction 1 — X(s) 
coming from “soft” processes, and assume that in nuclear collisions “hard” processes are 
incoherent and thus scale with the number of collisions, whereas “soft” processes scale 

m 


with the number of participants 


we arrive at the following simple parameterization 
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= [(1 - X(s)) (Np^rt) + X{s) {N,,a)] 


dr] 


dr] 


( 5 ) 


which describes the data quite well. In the framework of perturbative QGD approach, one 
has to assume that the coefficient X(s) is proportional to the mini-jet production cross 
section, and thus grows with energy reffecting the growth of the parton distributions at 
small Bjorken x, X(s) ~ [a:G(x)]^, with x l/^s. 
















Figure 3. Centrality dependence of the charged particle multiplicity near mid-rapidity in An 
collisions at ^/s = 20,130,and 200 GeV, from [ [l7||. 
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Therefore one expects 19] that the cen¬ 


trality dependence should become increas¬ 
ingly steep as the a/s increases (for the lat¬ 
est development, see however [11). This 
increase is not seen in Fig.^, which in the 
lower panel shows that the ratio of the dis¬ 
tributions at a/s = 200 GeV and ^/s = 130 
GeV is constant within error bars. The 
almost constant ratio appears consistent 
with the prediction [ ^,[ M based on the 



Figure 4: Pseudo-rapidity distributions of charged 
particles from Au -I- Au collisions at y/s = 200 GeV 
(open circles), from [|2^. Solid line is the prediction 
based on parton saturation [ ^ , and dashed line is 
the multi-phase transport model calculation 
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ideas of parton saturation, where the in¬ 
crease of multiplicity stems from the run¬ 
ning of the QGD coupling constant deter¬ 
mining the occupation number ~ 1/as of 
gluons in the classical field. 

A very important result presented at this Gonference is the centrality dependence of 
charged hadron multiplicity at RHIG injection energy of ^/s = 20 GeV [ 0, Fig.|. One 
can see that the shape of the centrality dependence changes very little over a large energy 
range, in which the perturbative minijet cross section grows by over an order of magnitude. 
The prediction of the saturation model [ l^j is seen to agree with the data; this indicates 
the possibility that parton saturation sets in in heavy ion collisions already at moderate 
energies. 
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2.1.3. Rapidity distributions 

Distributions of the produced particles in the emission angle 6 (with respect to the 
collision axis), or pseudo-rapidity rj = — ln[tan( 0 / 2 )] provide another important char¬ 
acteristic of the collision process. Two features of RHIC results (see Fig. are es¬ 
pecially noteworthy: i) the distributions do not exhibit scaling in rj] ii) the deviation 
from NN results is maximal in the central rapidity region whereas the shapes of the 
A A and NN distributions are similar in the fragmentation region. Moreover, when cor¬ 
rected for the different beam rapidity rj ^ rj — ybeami distributions at different ener¬ 
gies exhibit approximate scaling in the fragmentation region (“limiting fragmentation”). 


It will be very interesting to measure ra¬ 
pidity and centrality dependence of hadron 
production at RHIC in dA collisions. One 
of the really important questions which can 
be answered this way is whether there is 
a signihcant multiplication of the initially 
produced partons due to thermalization in 
Au — Au collisions [0], or the number of the 
produced hadrons is directly determined 
by the number of partons liberated early 
in the collision from the wave functions of 
the colliding nuclei. Very important con¬ 
straints in this respect are provided by the 
measurements of the total transverse en¬ 
ergy [ m . 



Figure 5: Azimuthal anisotropy of hadron produc¬ 
tion in Au -I- Au collisions at at ^/s = 130 GeV; V 2 
is the weight of the second harmonic, cos 2(p, in the 
particle distribution in the azimuthal angle (/?; from 

lE- 


2.1.4. Azimuthal distributions 

Of great interest and importance are the 

distributions of the produced hadrons in the azimuthal angle. Indeed, if all of the NN 
collisions were independent, there would be no reason to expect asymmetry in the dis¬ 
tribution of the produced hadrons in the azimuthal angle (measured with respect to the 
reaction plane). 

This is why the azimuthal asymmetry represents a sensitive test of the collective effects 
in nuclear collisions. The azimuthal angle distributions of the produced hadrons are 
usually expanded in harmonics in the following wayQ: 


dN 

dip 


1 -b 2vi cos if -f- 2 v 2 cos 2ip + ... 


( 6 ) 


Fig- I shows the extracted from RHIC data coefficient V 2 (v 2 7 ^ 0 in the language of the 
held corresponds to “elliptic how”). One can see that the asymmetry of the azimuthal 
distribution is quite sizable, and for peripheral collisions (small multiplicity Uchl'^max) 
reaches about 35%. 


^The absence of the terms proportional to sinn(/? is the consequence of parity conservation; it would be 
interesting to search for their presence in the data in view of the speculative scenarios allowing for P and 
CP violation in hot QCD [ |^ . 
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This effect certainly indicates the presence of collectivity in nuclear collisions, and comes 
out quite naturally in hydrodynamical calculations which assume complete thermalization 
in the hnal state [ ^], [ ^|- However, hnal state effects are not the only possible origin 


of the azimuthal asymmetry; indeed, as we have discussed above, the high degree of 
coherence in the initial state signaled by the multiplicity measurements, in the parton 
saturation scenario, introduces strong correlation between the transverse momentum of 
the parton and its transverse coordinate in the wave functions of the nuclei. When 
the nuclei collide, this effect mimics at least a part of the asymmetry usually ascribed 
exclusively to hnal-state interactions [ [0. The magnitude of the elliptic flow which 

can originate from the coherence in the initial state is still a subject of ongoing research 
and debates. 

The dependence of the elliptic flow on 
the transverse momentum of the hadrons 
presents a puzzle [ the value of V 2 is 


seen to hrst increase with pt, and then sat¬ 
urate. This contradicts to hydrodynam¬ 
ics, predicting a monotonic increase of V 2 
with Pt] of course hydrodynamics cannot 
be trusted at high pt anyway because the 
density of hard particles is too small to 
allow for a meaningful statistical descrip¬ 
tion. Energy loss of the produced jets can 
contribute to the azimuthal asymmetry at 
high Pt (see, e.g., [ 0), even though it is 
not yet clear if the magnitude of the effect 
can be reproduced under realistic assump¬ 
tions about the density of the medium and 
the jet interaction cross section [ ^ . 
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Figure 6: The yield of net protons at mid-rapidity 
in central collisions at different energies (black cir¬ 
cles); also shown is the inclusive yield of negative 
hadrons (open circles); from [ Q. 


2.1.5. Hadron abundances 

The measurements of yields of different hadrons at RHIC hold many surprises. Of 
particular interest is the fact that even at RHIC energies the asymmetry between baryons 
and antibaryons is still sizable, with p/p ratio about 0.65 [ |S3 


see Fig.^ This signals the 
- 10 - 2 . 


diffusion of baryon number to quite small values of Bjorken x 

Thermal statistical models are remarkably successful in explaining the yields of different 
hadrons at RHIC (see talks at this Conference, [|^, ^). This is certainly consistent with 
what is expected if the thermalized matter is produced. Can this success be considered 
a proof of thermalization? In my opinion, before we can conclude this we would have to 
understand better the origin of the successes of this model in explaining the particle yields 
in elementary e+e” and pp collisions, where conventional thermalization mechanisms are 
unlikely to operate. 


2.2. Hard processes 

2.2.1. Suppression of high pt particles 
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Jet energy loss was among the first sig¬ 
natures suggested for the diagnostics of the 
hot quark-gluon matter [ ^ . This 

is why the measurements of the high pt 
hadron production excited a lot of inter¬ 
est. Indeed, the experimental results are 
striking - as can be seen in Fig.^ the yield 
of high pt hadrons is drastically reduced 
with respect to what is expected for inco¬ 
herent production in NN collisions. This 
behavior is very different from what was 
previously observed in Pb — Pb collisions 
at CERN SPS and in a — a collisions at 
CERN ISR (see Fig.0). Does this impor¬ 
tant discovery signal jet energy loss in the 
quark-gluon plasma? The answer to this 
question can be given after we know the re¬ 
sults of the forthcoming measurements in 
p{d)A collisions, which would allow to sep¬ 
arate clearly the effects coming from the 
initial state. 
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Figure 7: The ratio of transverse momentum dis¬ 
tributions of charged hadrons and neutral pions in 
Au -I- Au and pp collisions at ^/s = 130 GeV; from [ 


2.2.2. Azimuthal correlations 

A very interesting recent result obtained at RHIC is the gradual disappearance of the 
back-to-back azimuthal correlations of high pt particles with centrality of the collision 
[ [ 1^, see Fig.||. This effect can be explained by the absorption of one of 


the jets in hot matter produced in central collisions. Another possible explanation stems 
from the onset of parton saturation in central collisions; this favors 2 —> 1 parton fusion 
processes as the leading production mechanism in central collisions, and such processes 
produce uncorrelated mono-jets. Again, a dehnitive answer regarding the origin of this 
intriguing effect will be given by the measurements in d—A collisions - while the saturation 
mechanism would still operate there, the absorption mechanism would not. 


2.2.3. R/vr puzzle 

Another striking puzzle at RHIC is a rapid increase of the baryon-to-pion ratio in 
central Au — Au collisions at high pt [ see Fig.^. The growth of this ratio is expected 
in the hydrodynamical scenario, in which equal velocity of the expanding parton “fluid” 
implies higher transverse momentum for more massive particles. However, the validity 
of hydrodynamical description is dubious at high pt where the density of particles is too 
small. If we assume that minijet fragmentation is the leading production mechanism of 
high Pt particles, then the growth of the R/tt ratio implies that minijet fragmentation is 
severely affected by the medium. 

Another scenario [ attempts to explain both R/vr puzzle and a large value of baryon 
asymmetry R/R ^ 1, and invokes the contribution of non-perturbative gluonic junctions 
in nuclear collisions 
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Figure 8. Azimuthal correlations of charged particles relative to a high pt trigger particle for peripheral 
(left) and central (right) collisions; from [Q, [ |4^. 



Figure 9. The ratio of antiprotons to pions as a function of transverse momentum for central and 
peripheral Au — Au collisions at 200 GeV (left), and the ratios of identified charged hadrons to neutral 
pions (right); from [ 


2.2.4. Charm production 

The production of heavy flavors and quarkonia represents a very important and exciting 
part of RHIC program. While these studies will benefit in the future from increased 
luminosity and improvements in the detectors (allowing, in particular, to reconstruct the 
decay vertex of heavy hadrons), the hrst measurement of charm production cross section 
has been already reported [ @]. This has been done by deciphering the charm decay 
contribution to the single electron spectrum - see Fig.p)0|. 

Of particular interest is the fact that while the production cross sections of light hadrons, 
as discussed above, show strong nuclear effects, the cross section of charm production, 
within the error bars of the measurement, scales with the number of iVAf collisions. These 
results may imply much smaller, in comparison to light quarks, energy loss of heavy quarks 
[ @1 stemming from the suppression of the gluon radiation at small angles (“dead cone 


effect”). 

Very intriguing hrst results on the production of charmonium have been reported at 
this Conference [ Since charmonium remains one of the best tools for the study 


of deconhnement [ |^, it will be extremely important to establish the dependence of 
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charmonium production cross section in nuclear collisions on centrality (for a discussion 
of various theoretical approaches to charmonium suppression at SPS, see e.g. [ Q). 


3. What have we learned so far? 
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It is too early to assess the implications 
of RHIC results; however, it is becoming 
increasingly clear that they challenge most, 
if not all, of the existing theoretical dog¬ 
mas. A coherent and convincing theory de¬ 
scribing all of the observed phenomena di¬ 
rectly in terms of QCD still has to be born. 

However, we can already conclude that 
many of the observed phenomena clearly 
manifest collective behavior; nuclear colli¬ 
sions at RHIC are not an incoherent super¬ 
position of nucleon-nucleon collisions. 

The measured particle multiplicities and 
transverse momentum spectra allow to es¬ 
timate initial energy density at the early 
moments of the collision; a typical value 
inferred in this way is about 20 GeV/fm^ 

(see, e.g., [ |^). The dynamics of strongly interacting matter at such energy density 
(exceeding the energy density in a nucleus by over two orders of magnitude!) should be 
described in terms of quarks and gluons, and the collective phenomena observed at RHIC 
thus directly reflect the properties of high density QCD. 


2 5 3 

‘ P.^(GeV/c) 


Figure 10: The background-subtracted electron 
spectra for minimum bias (0 — 92%) and central 
(0—10%) collisions compared with the expected con¬ 
tributions from open charm decays; from [ . 
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